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1 Introduction 

In 1,2 ' 3 , and 4 we presented a local mathematical model for EPR-type exper- 
iments. Our model is in agreement with the results predicted by Quantum 
Mechanics which in turn were confirmed by large scale experiments, first by 
Aspect, Roger and Dalibard 5 and later by several other teams. Due to space 
limitations the presentation of our model was rather terse in places. The 
purpose of this paper is to present our model in much greater detail and 
at the same time mathematically simplified because concerns and questions 
have been raised about non-locality and parameter dependence in a few recent 
publications 6 < 7 > 8 < 9 . 10 . Although we have answered these concerns at various 
occasions 11,12,13 , a comprehensive detailed exposition might be a better way 
to address these concerns. 

At first we give a brief summary of our model. In EPR type experiments 
two particles having their spin in a singlet state are emitted from a source and 
are sent to spin analyzers at two spatially separated stations S\ and S2. We 
assume with Bell that the particles emitted from the source are permitted to 
carry information in form of arbitrary hidden parameter random variables A 
that can assume values in some abstract space. 

In the original experiment by Aspect 5 , pairs of photons were emitted 
from the source once every few microseconds over a hour period governed by 
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a random process. We model this process mathematically in the following 
way. Imagine the time axis wrapped around a circle of circumference that 
corresponds to a time interval related to a simple measurement and normalized 
to 1. We suppose that for a fixed N each interval [(to — 1)/N,m/N], m = 
1, 2, . . . , N of arc length l/N on the circle gets about its proper share of time 
measurement points over the measurement period. This assumption, in turn, 
induces a random variable R, that we call the labelling variable, which assumes 
the values to = 1,2, ... ,N with equal probability, i.e., 

P(R = m) = l/N to = 1,2,..., N. (1) 

In Section 4 below, we present better motivation for the generation of the labels 
to and the random variable R by means of the Poisson process, commonly used 
to model spontaneous emissions. 

In our papers 1,2,3 ' 4 , using notation standard in Bell type proofs, the source 
parameter now denoted by A was not assigned a separate letter. The letter 
A, which in the standard notation generically symbolizes randomness, was 
used instead. As a consequence the label to being a function of the random 
emission times was misinterpreted by various authors as a function of A. In 
the present paper we have therefore decided to use the standard probability 
notation instead: Random variables are given separate names, they will be 
denoted by capital letters, they are measurable functions of uj or A, attached 
to some experiments, and there will be a clear distinction between the random 
variables, the values they can assume, and the set of measured data. 

After a pair of particles has been emitted from the source, the time of 
emission and thus of measurement are known, and so is the interval [(to — 
1)/N,m/N] on the unit circle into which the time of measurement falls. This 
determines the label to. While the pair of particles travels to their designated 
analyzer stations the experimenters (or a random number generator) can ex- 
ercise their free will and choose in their respective stations the directional set- 
tings, say a in Si and b in S2. Our model calls for hidden parameter random 
variables 1 A* t in S\ and A^* t in S2 which depend on the respective settings 
and on the time of measurement. The time of measurement is known (either 
the same or connected by a linear relation) in both stations. Consequently, 
the label to is known at both stations. This provides for the time correlation 
we elaborated in papers 1,2 ' 3 and 4 . 

Thus we have four random variables in operation, R, A, A* t , and A^* t . The 
joint density of all these variables p a b, is permitted 15 to depend on the settings 
a and b, and is given in Eq. (29). As a consequence of our construction we 
obtain certain properties of stochastic dependence relations which we state as 
a preview. In the general case where the distribution of A may depend on 
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time we have the following stochastic dependence relations between these four 
random variables signifying their time correlations. 

(i) The random variables A* t and A*,*, are stochastically independent. 

(ii) Given the random variable R the pair (A* t , A*,*.) is conditionally indepen- 

dent of A = A t . 

As to the probability distributions our construction yields the following prop- 



(iii) The probability distribution of A t can be chosen arbitrary. 

(iv) The probability distributions of A* t and A*,*, do not depend on a, b, nor 
t. 

In the special case where the distribution of A does not depend on time we 
have in addition to the above properties 

(ii)* The random variables A, A* t , and Aj^*. are stochastically independent. 

(vi)* The random variables R and A are stochastically independent. 

The random variables A a = ±1 and By, = ±1 symbolize the possible spin 
values and are functions only of a, A, R, and A* t , and of b, A, i?, and A*,*., 
respectively, thus obeying Einstein locality; that is 



erties. 



A a = A a {A* at ,A;R) 



and 



B h = B h (At u A;R). 



Moreover, we have 



E{A a B h } 



a b = 



cos <)(a, b) 



(2) 



and thus we have with probability 1 



B, 



(3) 



In addition, we have with probability 1 



E{A a | A, A* J = E{B h | A,Aj^} = 0. 



(4) 
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Note that the integration is, in essence, only performed with respect to R. If, 
in addition we integrate Eq. (4) with respect to A* t and AJjj we obtain with 
probability 1 

E{A a | A} = E{B h | A} = 0. (5) 

Since A & = ±1 and B\> = ±1, Eq. (5) implies parameter independence. How- 
ever, our variables A a and B\> depend on time through their functional depen- 
dence on R, a feature that is not considered in the original Bell definition of 
A a and £?t>- 

Thus our model is distinguished from standard Bell-type models by the 
introduction of the time related labelling variable R. If in Eq. (5) we also 
condition on R, besides conditioning on A, these equations no longer will hold. 
Moreover, we note that in the general case we have: 

(v) Given the random variable R the random variables A* t and A^ are 
stochastically dependent. 

(vi) The random variables R and At are stochastically dependent. 

(vii) The conditional probability distributions of A* t given R and of A£ t given 
R depend on both settings a and b. 

On a more basic level probabilities, conditional probabilities or even condi- 
tional expectations, such as the one in Eq. (5), can be interpreted as long term 
averages of outcomes of certain experiments. These long term averages can be 
thought of being taken over certain points on the time axis. We separate this 
averaging process into two parts by introducing the random variable R. We 
first average over the concatenated time intervals associated with a fixed label 
m = 1,2, ... ,N. Subsequently, we average the first averages over the values 
m that R can assume to obtain the overall averages. 

As a consequence we do not view the conditional stochastic dependence in 
(v) , nor the dependence on the settings of the conditional probability distribu- 
tions in (vii) as a violation of Einstein locality. These dependencies only signify 
the time correlations between the events in stations Si and S^- To express this 
in physical terms we point the reader to the following facts. The label m 
represents a concatenation of short time segments and not a given time, m 
therefore does not relate to or permit any instantaneous signalling. It can not 
be influenced by the experimenter in any significant way since it depends on 
the random spontaneous emission times and the largely arbitrary way of con- 
catenating these short time intervals in a specific interval [(m — 1)/N,m/N]. 
Therefore m is not an element of reality as opposed to, for example, the source 
parameter A. 
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Any argument for instantaneous action at a distance involving probabili- 
ties conditional on {R = m} must therefore be counterfactual. As an impor- 
tant example, one could argue that instead of a setting pair a, b the experi- 
menters might have chosen a, c. Then, since the joint probability conditional 
on {R — m} depends on both settings, the marginal distribution of A* t for 
setting a conditional on {R = m} may be different. How can this be without 
instantaneous action at a distance? The answer is that if c would have been 
chosen, then over a whole sequence of measuring times all the settings would 
be different. In order to have setting b with equal probability to setting c, 
the experimenters would have had to decide to choose c instead of b at other 
occasions. In other words, the whole history of settings would have to be dif- 
ferent. Because all the involved parameters, as well as the possible outcomes 
for the spin pair values may depend on the history, the probability distribution 
of A* t conditional on {R = m} may depend on the history and can therefore 
be different for the setting pairs a, b and a, c. The EPR argument postulates 
a physical reality of the source parameter A; in our papers we postulate also 
physical reality for the station parameters A* t and Aj^. However, we do not 
attach a physical reality in the same sense to the labelling random variable R. 
Fulfillment and violation of Einstein locality with respect to random variables 
such as the labelling variable R becomes a highly complex problem 16 . Let us 
note, in passing, that the exclusion of setting dependence conditional to any 
concatenation of time segments such as represented by m will automatically 
also exclude the result of the actual experiments which can be regarded as 
performed by concatenating the results obtained in certain time segments. 

We would like to emphasize that the joint probability measure given by 
Eq. (29) below is not canonical, i.e., not unique. This makes the model highly 
flexible to accommodate other possible set-ups of experiments. In fact, we hope 
to show with our work that the choice of the particular form of variables is 
mathematically highly flexible and can go far beyond simple ideas of elements 
of physical reality. We do not claim that the particular model actually exists in 
nature. All we want to show is that Bell type proofs actually can not do justice 
to the complexities involved in EPR experiments and therefore can not be 
used to draw conclusions about nonlocal effects as epitomized by instantaneous 
action at a distance. 

For clarity of presentation we develop our construction in several steps. 
The first one, almost identical with the presentation in ' , will be given in 
Section 2. In Section 3 we define the probability distribution of the time and 
setting dependent station parameters A* t and Aj^ . The measure we construct 
in these sections is not quite a probability measure (see Eq. (8), below). How- 
ever, it is a routine exercise to derive from it a probability measure by applying 
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some basic facts from the theory of weak convergence of probability measures. 
We may present the details in a paper to be submitted to a mathematics 
journal. 

We hope that this introduction provides enough of a guiding line through 
the mathematical intricacies that will follow. 



2 The First Step in Establishing the Model 

Before we start with the mathematics, let us recall that a pair of particles 
has been emitted from the source. The emission time and thus the measuring 
time is known. As a consequence the value m of the labelling variable R is 
determined. The experimenters have subsequently chosen their vectors a in 
S\ and b in S2, respectively. In effect, we assume that the measuring time, 
considered as a random variable and the labelling variable R are independent 
of the choice of vectors a and b. 

Let a = (ai, a 2 , a 3 ) and b = (61, 62, 63) be unit vectors. Our goal is to show 
that under our generalized conditions, it is possible to obtain the quantum 
result, the scalar product —a - b for the spin pair expectation value E{A CL B\ ] \. 
Here we formulate a theorem which provides the first stepping stone for this 
procedure. 

We define functions A a and B\> and choose the underlying measure space 
(R 2 7 1Z 2 ), i.e., the Euclidean plane {(u, v), —00 < u, v < 00} with Borel mea- 
surability, symbolized by 1Z 2 . We set 

!sign(a fe ) if - k < u < -k + 1 k = 1, 2, 3 
-1 ifj<w<j + i .7=0,1,... f . 

+ 1 if j + i <«< j + 1 .7=0,1,... (6) 

+ 1 elsewhere. 

Thus A depends here on a and u only. We will return below to the complete list 
of dependencies which only here would complicate the notation and not add 
to the present purpose. Here and throughout, we set sign(0) = 1. Similarly, 
we define 

-sign(& fe ) if - k <v < -k + 1 k = 1,2, 3 
„, , +1 Xj<v<j + \ .7=0,1,... 

bij 1-1 ifj + ±<v<j + lj = 0,l,... [7) 

-1 elsewhere. 

As in the case for A above B depends for the moment on b and v only. We 
now formulate the first step as a theorem. 
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Theorem 2.1 Let n > 4 be an integer. Then there exists a finite measure 
/i = /Ug'b with the following properties: jj, depends only on n, a and h, has 
compact support fl, satisfies 

1 < p{R 2 ) < 1 + 1/n 2 (8) 

(n) 

and has a density p = /cr b with respect to Lebesgue measure. Further 

[ A a (u)B h (v)p { ^(u,v)dudv = -ab (9) 

and for each vector a the following equation holds for all x: 

B a {x) = -A a {x). (10) 

The proof of the theorem requires the following fact which follows from a 
basic theorem on _B-splines 17 . We state the fact here in form of a lemma. 
Lemma 2.2 Let n > 4 be an integer. Then there exist real-valued functions 
Ni(x), ipi(y) with 1 < i < n depending only on real variables x and y, respec- 
tively, such that 

0<Ni(x)<l, 0<ipi(y)<2for0<x,y<l (11) 

and 

n 

< My)N t (x) -(y- xf < -n- 2 for0<x,y< 1. (12) 

i=l 

The proof of this lemma is given in Appendix 1. We now proceed to prove 
Theorem 2.1. 



Proof of Theorem 2.1: We first observe that Eq. (10) follows from the 
above definitions of A a and B^. Let il = [—3, 'in) 2 and let k be the indicator 
function of the union of the unit squares Uf"_ 2 [z — l,i) 2 , lined up along the 
main diagonal of ft, in symbols 

k(u, v)= ^ Sij - l{i-l<u<i}-l{j-l<v< j}. (13) 

i,j=-2, — l,...,3n 

Here 1{-} denotes the indicator function of the set in curly brackets and 

^-{SUiJi (14) 



7 



denotes the Kronecker symbol. On each of these in + 3 unit squares we place 
uniform mass, that comes from a product measure on each of the squares, 
where the first factor only depends on the setting a and the second factor only 
depends on the setting b. Although this will make the mathematics quite a bit 
more complicated, we can envision further experiments where this feature of 
our construction may be of importance. The details are as follows. We define 

\cik\ ■ l{-k < u < -k + 1} k e h 
N k (\ai\) ■ l{k- 1 < u < k} keh 
v*{u) = N k _ n (\a 2 \) -l{k-l<u<k} kel 3 (15) 
N k -2n(\a 3 \) -l{k-l<u<k}keh 
elsewhere 



\h\ ■ H~k < v < -k + 1} k e h 
lM\h\) ■ l{k - 1 < v < k} keh 
T b (v)= ±4> k _ n (\b 2 \) ■ l{k - 1 < v < k} kel 3 (16) 
\^k-2n{\h\) -l{k-l<u<k}keh 
elsewhere. 

The symbols Ji, . . . , I4 stand for I\ = +3, +2, +1; h = 1, . . . , n; ^3 = n + 
1, . . . , 2n; I4 = 2n + 1, . . . , 3n. We finally define the density by 



P { ab( u ' v ) = <y a (u)T b (v) K (u,v) (17) 

and the measure p by having density p^ with respect to Lebesgue measure. 

Hence we obtain from the above definitions the following integrals needed 
for the calculation of the spin pair correlation function: 



L 



A a (u)B b (v)p SLh (u,v)dudv = 

3,0) 2 

3 

|a fe ||6 fe |sign(a fe )sign(6 fe ) = a • b. (18) 

fe=i 



Furthermore, the integral over the complement of the square [—3, 0) 2 vanishes, 



i.e., 



/ A a (u)B\ 3 (v)p a b(u,v)dudv = (19) 

Jn\[-3,o) 2 



which proves Eq. (9). 

It remains to be shown that p a b defines a measure \i that is close to 
a probability measure, i.e., fulfills Eq. (8). For this, we consider the mass 
distribution between the square [—3, 0) 2 and its complement. The amount of 
mass Mi distributed over [—3, 0) 2 is 

3 

M 1 = Y,\"k\\bk\. (20) 
k=i 

The mass M 2 of fi\[-3,0) 2 equals 

-. 3 n 
fc=l i=l 

Thus the total mass distributed equals in view of Eq. (12) 

3 3 n 

M 1 +M 2 = ^ |a fc ||6 fc | + - £X>i(|a fc |)Vi(N) 
fe=i fe=i j=i 

3 3 

Mr + M 2 = MM + ^ Ed flfc l - N) 2 + 9 ■ n ~ 2 

fe=i fe=i 

Mi+M 2 = l + (9-n~ 2 (22) 

where < 6 < 1/4. 

This completes the proof of the theorem which is the first stepping stone 
of our construction of a suitable probability measure. 

Obviously, if instead of Eq. (21), we would define 

3 

M 2 := ^(|a fe |-|6 fc |) 2 , 

K=\ 

f2 = [—3, 3] 2 and would place the mass represented by these three summands on 
any of the nine unit squares of [0, 3) 2 , we would produce a genuine probability 
measure, satisfying all conclusions of Theorem 2.1. 

Finally, let L > 1, and pi > with J2e=i Vt = 1- F° r < w < 1, we define 

S (w) = (-iy, <=i 1 = 1,. 

q(w)=p t , f ^<w<{, £=l,...,L. (23) 
9 



For -3 < u, v < 3n and < w < 1 we define Q = ft x [0, 1) 



A a (u,w) = A a (u)s(w) 
B h (v,w) = B h (v)s(w) 
p ah (u,v,w) = p ab (u,v)q(w). 



Then A a and B b only depend on a, u, w, and b, v, w, respectively. Moreover, 
they satisfy the properly modified conclusion of Theorem 2.1. This procedure 
extends f2 to Q by adding as a factor the unit interval < w < 1 with a given 
mass distribution. 

3 Definition of the Layers 

We call the construction including the unit interval as factor, given in Section 2 
the first layer. To simplify the notation we shall omit the ~ sign from the fl, 
A, B, and p. As we noted in 1 and 3 the first layer does not yet provide 
a model that guarantees absence of action at a distance. To achieve this 
goal we will now define a system of layers. These layers will be obtained by 
permuting all the unit squares contained in fi, including the mass distribution 
and the corresponding strips on which A and B are defined. In addition we 
shall duplicate the mass distribution of each layer labeled m, labelling the 
duplicate layer m! . On the layer labeled m the functions A a m ^ and B h m ^ 
will remain unchanged. However, on the companion layer labeled m' we shall 



switch the signs of A a and B bl by setting Aa™ 1 ' = — A a m ^ and B h m •* — —B^ n \ 



As we observed in a recent paper 13 , this simple modification of our original 
construction encompasses all the desired features to achieve so called parameter 
independence. We now present this program in detail. 

Think of each of the unit cubes [i— 1, i) x \j— 1) x [0, l),i, j = —2, —1, . . . , 3n 
together with their respective mass distribution and the values of A a and By, 
defined on them as a unit ensemble. We permute these unit ensembles in 
the following way. Choose three vertical strips [i — x [— 3, 3n) x [0,1), 
i = — 2, — l,...,3n and three horizontal strips [— 3, 3n) x [j — l,j) x [0,1), 
j = —2, —1, . . . , 3n. These intersect in nine unit cubes. Place the three unit 
ensembles [i — x [i — x [0,1), i = —2,-1,0 of the first layer onto 
three of these nine unit cubes, such that each vertical and each horizontal 
strip contains exactly one of these three unit ensembles of the first layer, and 
move with them the vertical and horizontal strips of the first layer. This can 



assigned their unit ensembles. Choose 3n of them and place on them these unit 
ensembles of the first layer where the density was defined by ^Ni(\a,k\)ipi(\bk\), 




There are still 9n 2 unit cubes left to be 
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i = 1,2, ... , n; k = 1,2,3. This can be done in ^ 9 3 ™ (3n)! different ways. 

Place the remaining 9n 2 — 3n unit ensembles to fill up the empty spaces. They 
have total mass 0. This yields a grand total of 

arrangements, which we call "layers" . We call this number ^N. 

At this point we exercise our option to let — p m £ in Eq. (23), m = 
1,2,..., 7}N depend on the label of the layer. 

In summary, on each layer the functions A a (u, w; m) and B^v, w; m) only 
depend on (u,w) and {v,w), respectively. Each layer supports a measure 
Mm = Mabm satisfying 



l<Mm(lR 3 ) <l + l/n 2 . 

Each measure p m has a density p(u, v, w; m) with respect to Lebesgue measure 
that can be written in the form 

p(u, v, w; m) — <r a (u; m)rb(w; m)n(u, v; m)q(w; m) 

for —3 < u, v < in, < w < 1, with the obvious interpretation of <r a , Tb, k 
and q. Moreover, by Eq. (9), we have for each m = 1, . . . , \N , 

I A SL (u,w;m)B\ 3 (v,w;m)p(u,v,w;m)dudvdw — — a • b. (25) 

As indicated at the beginning of this section, we shall duplicate each layer so 
that at the end we will have a total of TV layers. We renumber the original layers 
to by the odd positive integers, 2m — 1, say, m = 1, 2, . . . , ^N. The companion 
layer to the layer 2m — 1 will be assigned label 2to, to — 1,2,..., ^N. Each 
layer 2m — 1 and companion layer 2m will be assigned density previously 
denoted p(u,v,w;m). Each layer 2to — 1 will carry the functions, originally 
denoted A a (u,w;m) and Bb(v,w;m), whereas the companion layer 2m will 
carry — A^u^w^m) and — B^v , w; m) , instead. Thus after renumbering the 
functions A a and B^ and the densities p accordingly we have for all u, w, w, 
and all to = 1, 2, ... , 

A a (u, w; 2m — 1) + A a (u, w; 2m) = 

B h (v,w;2m- 1) + B h (v,w;2m) = (26) 

and 

p(u, v, w; 2m — 1) = p(u, v, w; 2m). (27) 
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Moreover, Eq. (25) continues to hold for all m = 1, 2, . . . , N. Of course, the 
equivalent effect had been achieved by adding a fourth dimension t and multi- 
plying the original functions A & and £?b by a Rademacher function r(t). This 
was done in Section 5.3 of our paper 1 . 

With all the mathematical objects properly in place we now finalize the 
second step of the construction of our model. The emission time of the z-th 
particle determines the measurement time and thus the label m where m — 
1, 2, . . . , N. Recall that the labelling variable R has uniform distribution over 
the integers m — 1,2, ... ,N, given by Eq. (1). 

Apart from the random variable R the construction so far is plain calculus 
in ffi 3 . Only now we do define a realization of the random variables A* t , 
A£j, and A t by defining the conditional density of A* t , A^, and A t given the 
random variable R by 

Prob(A; t G [u, u + Am), A** G [v, v + Av),A t G [w, w + Aw) \R = m) * 
= p(u,v,w;m)AuAvAw 
= cr a (u; m)Tb(w; rn)K(u, v, m)q(w; m)AuAvAw 
— oo < u, v < oo, < w < 1; m = 1, 2, . . . , N. 

This is the same as saying that the joint density of the four random variables 
A* t , A**, A t , and R is given by 

p(u, v,w;m) • — , — oo < u, v < oo, < w < 1, m = 1, . . . , N. (29) 

A few remarks are in order. First, in previous write-ups we have included 
mappings / and g, to accommodate more general random variables A* t and 
A b j. Obviously, this can be done here, too. Second, we changed the model by 
defining p a b to be the joint conditional density of A* t and A^ given R, rather 
than by defining p a b given by Eq. (17) to be the joint conditional density of 
the mixed parameters A at and A bt , given R, as was done in b 3 . This makes 
for a more streamlined presentation when the source parameter A t is taken 
into account since obviously A at and A bt , are functions of A t and the station 
parameters A* t and A b j, respectively, and thus cannot be independent of At 
(compare to condition (ii*) in Section 1). Hence, the expression for the joint 
density corresponding to Eq. (29) would be more complicated. 

We now discuss the stochastic dependence relations between the four ran- 
dom variables A at , A b ^, A t , and R that are direct consequences of Eq. (29). 
First, the joint density of A at and A^* t is given by 

1 N f 1 

— ^ <T a (u;m)T h (v;m)K(u,v;m) / q(w;m)dw. 

m=l ^° 
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(3n + 3) 2 I — ' ' 2 

v y \fe=i fe=i i=i 



Since the last integral equals 1, this reduces to, in view of Eq. (22), 

/ 3 3 3 \ 

Ei^n^i + ^EE^(i« fe |)^(i^i)J = 

Mi + M 2 _ 1 + 6 ■ n- 2 
(3n + 3) 2 ~~ (3n + 3) 2 

with < 6 < j. We conclude that the joint density of A* t and A b *. is ap- 
proximately uniform over the square [—3, 3n) and, as a consequence, equals 
the product of its two marginal densities which are themselves approximately 
uniform over the interval [—3, 3n). 
We conclude that approximately: 

(i) A* t and A b *. are independent random variables, and 

(iv) the distributions of A* t and Aj^ do not depend on a, b, and t. 

Moreover, summation over I = 1,2, ... ,L yields 

Prob(A at G [u,u + Au),(A£* G [v, v + Av)\R = to) 

= <J a (u: m)Tb(v; to)k(u, v. m)AuAv, m = 1, . . . , N, — oo < u, v < oo 

and integration over — oo < u, v < oo yields 

Prob(A 4 G [w, w + Aw) | R = m) ~ q(w; m)Aw < w < 1, m = 1, . . . , N. 
Thus by Eq. (28), we have approximately 

(ii) Given R the pair (A* t , A^) is conditionally independent of A t . 
Also, approximately, 

Prob(A* t G [u, u + Am), (A£* G [v, v + Act) | A t G [w, w + Aw),R = m) 

cr a (u; TO)rb(cr; to)k(u, v. m)j^q{w; m)AuAvAw 
Prob(A t e[w,w + Aw) \ R = m)P(R = to) 
= cr a (u; TO)rb(v; m)n{u 1 v; m)AuAv 

= Prob(A* t G [u, u + Au), (A** G [v, v + Av) \ R = m). 

(ii*) Further ii p m £ = Pi independent of m, i.e., if At and R are independent, 
then A* t , A b *., and At are independent random variables. 
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Moreover, we obtain for the pair correlation integral 

E{A a B h } := E{A a (A: t ,A t ,R)B h (At* u A t7 R)} 
1 N 

m— 1 

by Eq. (25) and Eq. (29). Since by construction (see Eq. (28)) 

A (2m-1) + A (2m) = ^ B (2m-1) + fl (2 m ) = ^ m=1)2 ,...,^, 

we obtain parameter independence first summing over m to obtain Eq. (4) 
and then by keeping the desired variables fixed and by integrating over the 
remaining ones. 

4 A Model Based on the Poisson Process 

The original experiment of Aspect et al. 5 took hours to complete. Currently, 
improvements of the technique have been accomplished by various teams of 
experimenters 18 and the length of time it takes to perform these experiments 
has been reduced substantially. The time between subsequent measurements 
is still limited by the recovery (essentially a random process) of the detectors 
between two measurements. 

From the logistical angle the present section is designed to replace the 
third paragraph of Section 1 and the parts of Sections 2 and 3 corresponding 
to it. Thus, overall, the present section is a variant of that part of the model 
dealing with generating the labels m. This will be done by considering the 
waiting times between consecutive "jumps" of a Poisson process. Since we 
are entering more advanced mathematical territory we present some of the 
relevant definitions and theorems in a basic form rather than to send the reader 
searching through the literature. 

We first recall a few definitions from the theory of uniform distribution 
modi. For more details see 19 and 20 . For a real number x, denote by [x] the 
integer part and by {x} = x — [x] the fractional part of x. Let ({xi))°^ 1 be a 
sequence of real numbers. For k > 1 and < a < [3 < 1 denote by 

A k (a,0) :=J2 1 ( a <{x l }<P) 

i<k 

the number of elements Xi, i < k such that their fractional part {xi} is con- 
tained in a given interval [a, (3) C [0, 1). The sequence ((xj)) is called uniformly 
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distributed modi if its discrepancy 



D k := sup 

a, 



as k — > oo, that is, if in the long run, each interval [a, [3) uniformly contains 
its proper share of points {xi}. Equivalently, we could define A k by wrapping 
the real axis around a circle of circumference 1 and count the number of hits 
a given interval [a, [3), now located on the circle, receives from the sequence 
((xi)), i < k, itself. 

The standard mathematical model for spontaneous emissions of particles, 
such as photons or electrons, is a Poisson process with intensity 1/6, say. 
The waiting times Ti between successive emissions are independent identically 
distributed random variables having exponential distribution with parameter 
O. The following theorem is a special case of Theorem 2 of H. Robbins 21 . 

Theorem 4.1 Let ((T,)) be the sequence of waiting times between consecutive 
jumps of a Poisson process. Then with probability 1 the sequence ((Ti + • • • + 
Ti))iZi is uniformly distributed modi. 



Remark In fact, it follows easily from Robbins' proof and the Erdos-Turan 
inequality that with probability 1 the discrepancy D k tends to zero at least 
with speed k~ 2 (log k) 3 . There are more than a dozen other papers extending 
Robbins' theorem. 

In terms of weak convergence of probability measures Theorem 4.1 can 
be reformulated in the following way. (See e.g., Billingsley 22 , pp. 15-25.) 
Let uj be an element of the set f2* of probability 1 as in Theorem 4.1. Set 

Xi = Ti (ui) H h Ti(u) . Let P k be the probability measure that assigns point 

mass j to each {xi} 7 1 < i < k. If several {x^ coincide, let the mass add. 
Then P k P in the sense of weak convergence. Here P denotes Lebesgue 
measure on [0,1). For ease of presentation let us reformulate Theorem 4.1 
in terms of random variables. Let Y k be a random variable defined on some 
probability space such that 

Prob(F fc = { Xi }) = j i = l,2,...,fc 
k 

and let U be a random variable having uniform distribution on [0, 1), i.e., 
Prob(C7 < x) = x, < x < 1. Then Theorem 4.1 can be restated as follows. 
For each w e ft* , we have as k — > 00 



Y k =>U 
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in the sense of weak convergence. Let N be defined in Eq. (24). For to — 
1,2, ... ,N define intervals I m , of length 1/N by 

m — 1 m\ 

m := [^T'nJ ■ 

Define a new random variable R by setting 

R = m if U e I m ,m= 1,2, ...,N. 
Then uniformly over all intervals I m , m = 1, 2, . . . , N, we have 

Prob(i? = m) = Prob([7 G 7 m ) = lim Prob(T fe G J m ) = (30) 

As in Section 3, suppose that the i-th pair of particles has been emitted. 
Fix w G 0*. The time of emission of the i-th pair equals a;, = Ti («;)+• • -+Ti{uj). 
When reduced modi the fractional part {xi} determines a label m with 1 < 
m < N . The labelling variable R has uniform distribution over the integers 
to = 1, 2, . . . , TV given by Eq. (30) or Eq (1). However, at the time the pair of 
particles arrive at their respective measuring stations, the devices may not yet 
be ready to provide a measurement, because of recovery problems, etc. For 
r = 1, 2 we define the random variable D r by setting 

D r = +1 if device is ready at station S r 
= if not. 

Obviously, by stochastic independence, 

P(R = m\D u D 2 ) = P(R = to) = m = 1, 2, . . . , N. 

Hence, given that the devices at both stations are ready for measurement, 
the labelling random variable R still has uniform distribution. 
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Appendix 1 

The lemma is an immediate consequence of Theorem 4.1 of Schumaker 17 for 
the special values of m = 3, i = 0, r = n, and the knots chosen to be y v = ^ 
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with v = 0, ±1, ±2, . . . Then by Schumaker's equation (4.33) we have (dropping 
the fixed superscript 3 of Nf): 

n 

(y - xf = </>i,3(v) N i( x ) for all < x < 1 and y G i?. 

i=-2 

Here 

4>iAv) = (y- Vi+i)(y - ^+2) 

and 

< Ni(x) < 1 for all x. 

We now restrict y to < y < 1. Then for — 2 < i < n, we have < <j>i,z(y) < 2 
unless y G y i+2 ]. Since we must avoid negative <f), we set 4> — in this 

interval by defining new functions tp: 

ipi{y) = if y £ [2/i+i,2/i+2] 
= 4>i^{y) otherwise. 

Since for y G [^+1,^+2] we have 

l(y-Wi+i)(y-yi+2)| < 

we have 

n , 

< J2 (MvWx) - 4>MNi{x)) < ^ 

i=-2 

because for any given y and for all x, only one term in the sum can be off by 
at most jij. This proves the lemma. 
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